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ABSTRACT 

Robust gear designs consider not only crack initiation, but 
crack propagation trajectories for a fail-safe design. In actual 
gear operation, the magnitude as well as the position of the 
force changes as the gear rotates through the mesh. A study to 
determine the effect of moving gear tooth load on crack 
propagation predictions was performed. Two-dimensional 
analysis of an involute spur gear and three-dimensional analysis 
of a spiral-bevel pinion gear using the finite element method 
and boundary element method were studied and compared to 
experiments. A modified theory for predicting gear crack 
propagation paths based on the criteria of Erdogan and Sih was 
investigated. Crack simulation based on calculated stress 
intensity factors and mixed mode crack angle prediction 
techniques using a simple static analysis in which the tooth load 
was located at the highest point of single tooth contact was 
validated. For three-dimensional analysis, however, the 
analysis was valid only as long as the crack did not approach 
the contact region on the tooth. 


INTRODUCTION 

Effective gear designs balance strength, durability, 
reliability, size, weight, and cost. Even effective designs, 
however, can have the possibility of gear cracks due to fatigue. 
In addition, truly robust designs consider not only crack 
initiation, but crack propagation trajectories. As an example, 
crack trajectories that propagate through the gear tooth are the 
preferred mode of failure compared to propagation through the 
gear rim. Rim failures will lead to catastrophic events and 
should be avoided. Analysis tools that predict crack 
propagation paths can be a valuable aid to the designer to 
prevent such catastrophic failures. 


Pertaining to crack analysis, linear elastic fracture 
mechanics applied to gear teeth has become increasingly 
popular. The stress intensity factors are the key parameters to 
estimate the characteristics of a crack. Analytical methods 
using weight function techniques to estimate gear tooth stress 
intensity factors have been developed (Nicoletto, 1993, 
Abersek and Flasker, 1994). Numerical techniques such as the 
boundary element method and finite element method have also 
been studied (Sfakiotakis, et al., 1997, Inoue and Kato, 1994). 
Based on stress intensity factors, fatigue crack growth and gear 
life predictions have been investigated (Blarasin, et al., 1997, 
Glodez, et al. 1998, Arikan, et al, 1998, and Abersek and 
Flasker, 1998). In addition, gear crack trajectory predictions 
have been addressed in a few studies (Lewicki and Ballarini, 
1997, Pehan, et al., 1997, Curtin, et al., 1998, Lewicki, et al., 
1999, and Ciavarella and Demelio, 1999). 

From publications on gear crack trajectory predictions, the 
analytical methods have been numerical (finite or boundary 
element method) while solving a static stress problem. In 
actual gear applications, however, the load moves along the 
tooth, changing in both magnitude and position. No work has 
been done investigating the effect of this moving load on crack 
trajectories. 

The objective of the current work is to study the effect of 
moving gear tooth load on crack propagation predictions. Two- 
dimensional analysis of an involute spur gear using the finite 
element method is discussed. Also, three-dimensional analysis 
of a spiral-bevel pinion gear using the boundary element 
method is discussed. A quasi-static numerical simulation 
method is presented in which the gear tooth engagement is 
broken down into multiple load steps and analyzed separately. 
Methods to analyze the steps are discussed, and predicted crack 
shapes are compared to experimental results. 


NASA/TM—2000-210227 


1 




TWO-DIMENSIONAL ANALYSIS 

Gear Modeling. The two-dimensional analysis was 
performed using the FRANC (FRacrure ANalysis Code) 
computer program developed by Wawrzynek (1991). The 
program is a general-purpose finite element code for the static 
analysis of tvvo-dimensional cracked structures. The program 
uses principles of linear elastic fracture mechanics and is 
capable of analyzing plane strain, plane stress, or axi- 
symmetric problems. A unique feature of the program is the 
ability to model cracks and crack propagation in a structure. A 
rosette of quarter-point, six-node, triangular elements is used 
around the crack tip to model the inverse square-root stress 
singularity. Mode I and mode II stress intensity factors, Kj and 
K/i, respectively, can be calculated using a variety of methods. 
(As a refresher, mode I loading refers to loads applied normal 
to the crack plane and tends to open the crack. Mode II refers 
to in-plane shear loading.) The stress intensity factors quantify 
the state of stress in the region near the crack tip. In the 
program, the stress intensity factors can be used to predict the 
crack propagation trajectory angles, again using a variety of 
methods. In addition, the program has a unique re-meshing 
scheme to allow automated processing of the crack simulation. 

A spur gear from a fatigue test apparatus was modeled to 
demonstrate the two-dimensional analysis. The modeled gear 
had 28 teeth, a 20° pressure angle, a module of 3.175 mm 
(diametral pitch of 8 in* 1 ), and a face width of 6.35 mm (0.25 
in). The gear had a backup ratio (defined as the rim thickness 
divided by the tooth height) of 3.3. The complete gear was 
modeled using mostly 8-node, plane stress, quadrilateral finite 
elements. For improved accuracy, the mesh was refined on one 
of the teeth in which a crack was inserted. The total model had 
2353 elements and 7295 nodes. Four hub nodes at the gear 
inner diameter were fixed to ground for boundary conditions. 
The material used w*as steel. 

Tooth Loading Scheme. To determine the effect of gear 
tooth moving load on crack propagation, the analysis was 
broken down into eighteen separate load cases (Tig. 1). An 
initial crack of 0.26 mm (0.010 in) in length was placed in the 
fillet of tooth 2, normal to the surface, at the location of the 
maximum tensile stress (uncracked condition). Six load cases 
were analyzed separately with the load on the tooth ahead of 
the cracked tooth, six on the cracked tooth, and six on the tooth 
after. The calculated stress intensity factors for unit loads at 
each of the load positions is shown in Fig. 2. These stress 
intensity factors were calculated using the J-integral technique 
(Rice, 1968). Loads on tooth 2 (cracked tooth) produced 
tension at the crack tip. The K/s increased as the load moved 
toward the tooth tip (load cases 12 to 7, Fig. 2b) due to the 
increased load lever arm. Loads on tooth 3 also produced 
tension at the crack tip, but at an order of magnitude less than 
those produced from the loads on tooth 2 (Fig. 2c). Loads on 
tooth 1 gave compression to the crack tip as shown by the 
negative K/s (Fig. 2a). 

Next, the actual load magnitudes on the gear tooth were 
considered as it went through the mesh. Computer program 
DANST (Dynamic ANalysis of Spur gear Transmission, Lin, et 


ah. 1988) was used for the analysis. This program is based on a 
four-degree-of-freedom, torsional, lumped mass model of a 
gear transmission. The model includes driving and driven 
gears, connecting shafts, a motor, and a load. The equations of 
motion for this model were derived from basic gear geometry, 
elementary vibration principles, and time-varying toorh 
stiffnesses. For simplicity, the static gear tooth loads of the 
solution were determined (Fig. 3). These loads were 
determined from well-established gear tooth stiffness principles 
and static equilibrium. The loads are shown as a function of 
gear rotation for a driver torque of 68 N-m (599 in-lb). Tooth 2 
began contact at a gear rotation of 10°. As the gear rotation 
increased, the load on tooth 2 gradually increased. Tooth 1 and 
2 shared the load fora rotation from 10 to 18°. From 18 to 25 
tooth 2 carried the complete load. At 23°, tooth 2 is considered 
at its highest point of single tooth contact (HPSTC). 

The stress intensity factors as a function of gear rotation 
were then determined by multiplying the stress intensity factors 
determined from the units loads (Fig. 2) by the actual tooch 
loads (Fig. 3) and applying superposition since linear elastic 
fracture mechanics was used. The results are showrn in Fig. 4. 
As expected, the mode I stress intensity factor (Fig. 4a) was 
mostly influenced by the load on tooth 2. Note that the largest 
value of Kj occurred at the HPSTC. Also note that the 
magnitude of Kj (Fig. 4a) w'as much larger than that of K I} (Fig. 
4b). This implied that Kj was the driving force in the crack 
propagation. K Ih however, affected the crack propagation angje 
as will be shown in the next section. 

Crack Propagation Simulation. From Williams (195" i. 
the tangential stress near a crack tip, is given by 

G( * = 7^( A ' cos' 7 - 3K * s^ f cos 2 1 j < 1 1 

where r and 6 are polar coordinates with the origin at the crack 
tip. Erdogan and Sih (1963) postulated that crack extension 
starts at the crack tip and grow's in direction of the greatesi 
tangential stress. The direction of the greatest tangential stress 
is determined by taking the derivative of Eq. 1 with respect to 
ft setting the expression equal to zero, and solving for ft 
Performing the math, this predicted crack propagation angle, 
ftm is given by 



From Eq. 2, the predicted crack propagation angle is a 
function of the ratio of K) to K //. Erdogan and Sih (1963) used 
brittle plexiglass plates under static loading to validate their 
proposed theorems (i.e., the ratio of K } to K u was constant). For 
the gear problem in the current study, however, the ratio of K : 
to K n was not constant during gear rotation. This is shown in 
Fig. 4c (actually plotted as the ratio K n to Kj for clarity). In 
addition. Fig. 4d gives the calculated 6 m from Eq. 2 as a 
function of gear rotation. 
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In order to simulate gear crack propagation, a modification 
to the Erdogan and Sih theory was postulated in the current 
study. This modified theory states that the crack extension 
starts at the crack tip and grows in the direction of the greatest 
tangential stress as seen during engagement of the gear teeth. 
The procedure to calculate the crack direction is as follows: 
1)A^/ and K u are determined as a function of gear rotation 
(Figs. 4a and 4b, as described in the previous section), 2) the 
ratio of Kj to K u as a function of gear rotation is determined 
(Fig. 4c), 3) 6 m (using Eq. 2) as a function of gear rotation is 
determined (Fig. 4d), 4) Oqq (using Eq. 1) as a function of gear 
rotation is determined (Fig. 4e), 5) the predicted crack direction 
is the value of G m for which Gee is greatest during gear rotation. 

For the gear example given, the tangential stress factor 
(defined as Gqq xr ) is plotted as a function of gear rotation 
in Fig. 4e. This plot looks very similar to the mode I stress 
intensity factor plot (Fig. 4a) since Kj was much larger than K u 
(see Eq. 1). The tangential stress was largest at the HPSTC 
(gear rotation of 23°) and the predicted crack propagation angle 
at this gear rotation was 6 m = 4.3°. 

Using this propagation angle, the crack was extended by 
0.26 mm (0.010 in), re-meshed, re-analyzed, and a new 
propagation angle was calculated using the method described 
above. This procedure was repeated a number of times to 
produce a total crack length of 2.38 mm (0.094 in). The 0.26- 
mm crack extension length was based on prior experience in 
order to produce a smooth crack path. Fig. 5 shows the stress 
intensity factors versus gear rotation for a number of crack 
lengths. Note that the mode I stress intensity factors looked 
similar but with increased magnitude as the crack length 
increased. In all cases, the selected crack propagation angles 
occurred when the tooth load was placed at the HPSTC. Fig. 6 
shows a similar analysis but with a model of a thin-rimmed 
gear. Here, the gear was modeled based on the previous design, 
but with slots incorporated in the rim to simulate a thin-rimmed 
gear. The backup ratio for this model was 0.2. As seen, the 
magnitudes of the mode I stress intensity factors during tension 
(gear rotations 18 to 45°) were larger than that of the 3.3 
backup ratio gear. Also, there was a significant increase in the 
compressive Kf s (gear rotation less than 18°) due to the 
increased compliance of the thin rim gear. 

Comparison to Experiments. Fig. 7 shows the results of 
the analysis compared to experimental tests in a gear fatigue 
apparatus. The original model (backup ratio of 3.3) as 
described before was compared along with models of backup 
ratio of 1.0 and 0.3. These later two models were created using 
slots in the gear blank as previously described. The 
experiments were first reported by Lewicki and Ballarini 
(1997). Here, notches were fabricated in the tooth fillet region 
to initiate tooth cracking of test gears of various rim 
thicknesses. The gears were run at 10,000 rpm and at a variety' 
of increasing loads until tooth or rim fracture occurred. As 
seen from the figure, good correlation of the predicted crack 
trajectories to experimental results was achieved. For backup 
ratios of 3.3 and 1.0, tooth fractures occurred. For the backup 


ratio of 0.3, rim fracture occurred. 

As a final note, the analysis indicated that the maximum 
tangential stress at the crack tip always occurred when the tooth 
load was positioned at the HPSTC. Thus, for two-dimensional 
analysis, crack simulation based on calculated stress intensitv 
factors and mixed mode crack angle prediction techniques car 
use a simple static analysis in which the tooth load is located a: 
the HPSTC. This was based on a modification to the Erdogar 
and Sih crack extension theory' and the fact that the mode I 
stress intensity factor was much larger than the mode II factor. 


THREE-DIMENSIONAL ANALYSIS 

Gear Modeling. The three-dimensional analysis was 
performed using the FRANC3D (FRacture ANalvsis Code for 5 
Dimensions) computer program developed by Wawrzvnek 
(1991). This program uses boundary element modeling anc 
principles of linear elastic fracture mechanics to analvze 
cracked structures. The geometry of three-dimensional 
structures with non-planar. arbitrary shaped cracks can be 
modeled. The modeling of a three-dimensional cracked 
structure is performed through a series of programs. Structure 
geometry grid point data are imported to a solid modeler 
program. Here, appropriate curves and faces (or patches) are 
created from the grid data as well as a closed-loop surface 
geometry model. This surface model is then imported to the 
FRANC3D program for boundary' element model preparation. 
The user can then mesh the geometry model using 3 or 6 node 
triangular surface elements, or 4 or 8 node quadrilateral 
elements. Boundary conditions (applied tractions anc 
prescribed displacements) are applied on the model geometry 
over faces, edges, or points. Initial cracks, such as elliptical or 
penny shaped, can be inserted in the structure. After complete 
formulation, the model is shipped to a boundary element 
equation solver program. Once the displacement and traction 
unknowns are solved, the results are exported back to the 
FRANC3D program for post processing. Fracture analysis, 
such as stress intensity factor calculations, can then be 
performed. 

The spiral-bevel pinion of the OH-58C helicopter main 
rotor transmission was modeled to demonstrate the three- 
dimensional analysis. The pinion had 19 teeth, a 20° pressure 
angle, a 30° mean spiral angle, a module of 3.66 mm (diametral 
pitch of 6.94 in 1 ), and a face width of 31."5 mm (1.28 in). For 
OH-58 operation, the pinion mates with a -tooth spiral-bevel 
gear, operates at 6060 rpm, and has a design torque of 350 N-m 
(3099 in-lb). 

The boundary element model of the OH-58 pinion 
developed by Spievak (1999) was used for the study. Three 
teeth, the rim cone, and the bearing support shafts were 
modeled (Fig. 8). The tooth surface and fillet coordinates were 
determined from the methods developed by Handschuh and 
Litvin (1991) and Litvin and Zhang (1991). The mesh of the 
three teeth was refined for improved accuracy. A half-ellipse 
initial crack with major and minor diameters of 3.175 and 2.540 
mm, respectively (0.125 and 0.100 in), was placed in the fillet 
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of the middle tooth normal to the surface. The crack was 
centered along the face width and centered along the fillet. The 
complete gear model had a total of about 2600 linear elements 
(both triangular and quadrilateral) and about 2240 nodes. For 
boundary conditions, the end nodes of the larger diameter shaft 
were fixed and the nodes on the outer diameter of the smaller 
diameter shaft were constrained in the radial directions. Again, 
the material was steel. 

Tooth Contact Analysis and Loading Scheme. Due to 
the geometrical complexities and three-dimensional action, 
numerical methods are required to determine the contact loads 
and positions on spiral-bevel teeth since no closed-form 
solution exists. The method of Litvin and Zhang (1991) was 
used to determine the mean contact points on the spiral-bevel 
pinion tooth. The method modeled tooth generation and tooth 
contact simulation of the pinion and gear. With the mean 
contact points taken as the centers, contact ellipses were 
determined using Hertzian theory (Handschuh and Richer, 
1996). Fig. 9 shows the estimated contact ellipses on the spiral- 
bevel pinion tooth. Fifteen separate ellipses (load cases) were 
determined, starting from the root of the pinion and moving 
toward the tooth tip and toe. Load cases 1-4 and 12-15 w f ere 
double tooth contact regions while load cases 5-11 were single 
tooth contact regions. Note that load case 11 corresponds to the 
load at the HPSTC. For each load case using the boundary 
element method, tractions were applied normal to the surface to 
the appropriate ellipse with the magnitude equal to the tooth 
normal force divided by the ellipse area. 

Crack Propagation Simulation, The procedure for the 
three-dimensional crack propagation simulation of the OH-58 
spiral-bevel pinion was as follows. For each of the load cases 
of Fig. 9, the mode I and mode II stress intensity factors were 
determined at 25 points along the crack front (note that for 
three-dimensions, there is a crack front, not just a crack tip as in 
two-dimensions). The extended crack directions at each of 
these 25 points were determined using the modified Erdogan 
and Sih crack extension theory as described in the two- 
dimensional analysis. That is, as the cracked spiral-bevel 
pinion tooth was engaged in the mesh, the crack extension 
started at each point along the crack front and grew in direction 
of the greatest tangential stress at those points during mesh. 
The amount of crack extension at each point along the crack 
front was determined based on the Paris crack growth 
relationship (Paris and Erdogen, 1963) where 


where a ,* was the amount of extension of the i th point along the 
crack front, K Iti was the mode I stress intensity factor of the i th 
point along the crack front corresponding to the load case 
which gave the largest tangential stress for that front point, 
£/.max was the value of the largest K u along the crack front, a ^ 
was the maximum defined crack extension along the crack 
front, and n was the Paris material exponent. From experience, 
the maximum extension size, was set to 1.27 mm (0.050 
in.). The Paris exponent, n, was set to 2.954 based on material 


tests for AISI 9310 steel by Au and Ke (1981). A third-order 
polynomial was then used to smooth the extended crack front. 
The new crack geometry was then re-meshed. After re¬ 
meshing, the model was rerun and solved for stress intensity 
factors, and crack propagation directions. The above procedure 
w as repeated a number of times to simulate crack growth in the 
gear tooth. 

Table I gives results from the first four steps during this 
process. Note that step 0 corresponds to the initial half-ellipse 
crack. For steps 0 and 2, the largest tangential stress occurred 
at the HPSTC (load case 11) for the majority of the points along 
the crack front. For steps 1 and 3, the largest tangential stress 
occurred at load cases 8, 9, 10, or 11. 

Fig. 10 shows the stress intensity factor distribution along 
the crack front for step 1 (crack area of 5.96 mm : (0.009 in 2 )). 
Similar to the spur gear analyses, the K/s were larger as the 
load moved from the root to the tip due to the larger load lever 
arm. Other than absolute magnitude, the Kj distributions along 
the crack front looked similar for the various load cases. Fig. 
11 depicts the stress intensity’ factors plotted against load case 
(at a point along the front, biased toward the toe, normalized 
position along the crack front of 0.83). This figure show's the 
simulated distribution as the pinion engages in mesh with the 
gear. Note again that the ratio of K n to Kj was not constant 
during engagement. 

Fig. 12 shows exploded views of the pinion crack 
simulation after seven steps. It should be noted that the loading 
was placed only at the HPSTC for the last three steps. This was 
due to modeling difficulties encountered using the multi-load 
analysis. It was felt that this simplification did not significantly 
affect the results due to the smoothing curve-fit used. In 
addition, the tangential stress near the crack tip was either 
largest, or near its largest value, when the load was placed at 
the HPSTC. 

Comparison to Experiments. Fig. 13 shows the results of 
the analysis compared to experimental tests. The experimental 
tests were performed in an actual helicopter transmission test 
facility. As was done with the gear fatigue tests described 
before, notches were fabricated in the fillet of the OH-58 pinion 
teeth to promote fatigue cracking. The pinion was run at full 
speed and a variety of increasing loads until failure occurred. 
Showm in the figure are three teeth that fractured from the 
pinion during the tests (Fig. 13b). Although the notches were 
slightly different in size, the fractured teeth had basically the 
same shape. 

A side view of the crack propagation simulation is shown 
in Fig. 13a for comparison to the photograph of the tested 
pinion in Fig. 13b. From the simulation, the crack immediately 
tapered up toward the tooth tip at the heel end. This trend 
matched that seen from the tests. At the toe end, the simulation 
show r ed the crack progressing in a relatively straight path. This 
also matched the trend from the tests. Toward the latter stages 
of the simulation, however, the crack tended to taper toward the 
tooth tip at the toe end. This did not match the tests. One 
problem encountered in the simulation during the later steps 
was that the crack at the heel end of the tooth became close to 


NASAtTM—2000-210227 


4 




the actual contact ellipses. It was felt that the crack-contact 
interaction may have influenced the trajectory predictions to 
cause the discrepancy. 

Spievak (1999) reported on another method to account for 
the non-uniform Kjj to A7 ratio during pinion tooth ensaeement. 
This method considered contributions from all load cases in the 
crack angle prediction scheme and presented a method to 
cumulate the load effects. From these studies, reported crack 
propagation simulation of an OH-58 pinion also predicted the 
erroneous taper toward the tooth tip at the toe end. Again, the 
crack-contact interaction may have influenced the trajectory 
predictions to cause the discrepancy. Spievak also reported on 
a simulation using only the load at the HPSTC. The results 
from that simulation were similar to that achieved in the current 
study. It should be noted that the proposed method in the 
current study to account for moving tooth loads for the three- 
dimensional analysis was extremely cumbersome. It is 
therefore felt that the analysis using only the load at the HPSTC 
appeared accurate as long as the crack did not approach the 
contact region on the tooth. 

CONCLUSIONS 

A study to determine the effect of moving gear tooth load 
on crack propagation predictions was performed. Two- 
dimensional analysis of an involute spur gear using the finite 
element method was investigated. Also, three-dimensional 
analysis of a spiral-bevel pinion gear using the boundary 
element method was discussed. The following conclusions 
were derived: 

1) A modified theory for predicting gear crack 
propagation paths based on the criteria of Erdogan and Sih was 
validated. This theory stated that as a cracked gear tooth was 
engaged in mesh, the crack extension started at the crack tip 
and grew in direction of the greatest tangential stress during 
mesh. 

2) For two-dimensional analysis, crack simulation based 
on calculated stress intensity factors and mixed mode crack 
angle prediction techniques can use a simple static analysis in 
which the tooth load is located at the highest point of single 
tooth contact. 

3) For three-dimensional analysis, crack simulation can 
also use a simple static analysis in which the tooth load is 
located at the highest point of single tooth contact as long as the 
crack does not approach the contact region on the tooth. 
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Step 

Crack area 
(mm 2 ) 

Crack front point(s) 

Load case for 
largest 

0 

3.12 

1 

8 



2-25 

ii 

1 

5.96 

1 

9 



2,4-7,21,23-25 

10 



3,8-20,22 

11 

2 

10.35 

1 

10 



2-25 

11 

3 

13.35 

7,9,20 

8 



5-7,10-15.21,26.2'' 

9 



1-4,8,16-19.22-25 

11 



Figure 1. Location of load cases for finite element mesh. 



0 5 10 15 20 25 30 35 

Tooth roll angle, deg 

Figure 2. Mode I and mode II stress intensity factors for a unit load 
and an initial crack of 0.26 mm. 



Figure 3. DANST computer program output of 
static gear tooth load, 68 N-m driver torque. 
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a) Backup ratio = 3.3. b) Backup ratio = 1.0. 



c) Backup ratio = 0.3. 

Figure 7. Comparison of predicted gear tooth crack 

propagation paths with experimental results 
(P = predicted. E = experiments). 


10 



Load case 


Figure 9. Location of tooth contact ellipses and magnitude of load 
on OH-58 spiral-bevel pinion tooth. 



Figure 8. Boundary element model of OH-58 spiral-bevel pinion 



Crack front 






Figure 10. Stress intensity factors from three-dimensional 
OH-58 pinion tooth crack propagation simulation; 
step 1, crack area = 5.96 mm 2 . 
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Normalized position 
along crack front = 0.83 



1 3 5 7 9 11 13 15 

Load case 

Figure 11. Stress intensity factors from three-dimensional 
OH-58 pinion tooth crack propagation 
simulation; step 1, normalized position along 
crack front = 0.83. 


Crack mouth 


Heel 



a) Simulation. b) Experiments. 

Figure 13. Comparison of OH-58 spiral-bevel pinion tooth 
crack propagation simulation to experiments. 




Figure 12. OH-58 spiral-bevel pinion tooth crack propagation 
simulation after seven steps. 


NAS A/TM—2000-210227 


9 






REPORT DOCUMENTATION PAGE 


Form Approved 
OMB No. 0704-0188 


Public reporting burden lor this collection of information is estimated to average 1 hour per response, including the time for rev>e*r>g instructions, searching existing data sources, 
gathering and maintaining the data needed, and completing and reviewing the collection of information. Send comments regan*rg r*s burden estimate or ary other aspect of this 
collection of information, including suggestions for reducing this burden, to Washington Headquarters Services. Directorate for Irfry-^on Operations and Repocs. 1215 Jefferson 
Davis Highway, Suite 1204. Arlington. VA 22202-4302. and to the Office of Management and Budget. Paperwork Reduction Poecr 0704-01861. Washington. EC 20503. 


1. AGENCY USE ONLY (Leave blank) 


TITLE AND SUBTITLE 


REPORT DATE 

Julv 2000 


REPORT TYPE AND DATES COVERED 

Technical Memorandum 


5. FUNDING NUMBERS 


Consideration of Moving Tooth Load in Gear Crack Propagation Predictions 


6. AUTHOR(S) 

David G. Lewicki, Lisa E. Spievak, Paul A. Wawrzynek, 
Anthony R. Ingraffea, and Robert F. Handschuh 


7. PERFORMING ORGANIZATION NAME(S) AND ADDRESS(ES) 

NASA Glenn Research Center 
Cleveland. Ohio 44135-3191 
and 

U.S. Army Research Laboratory 
Cleveland. Ohio 44135-3191 


9. SPONSORING/MONITORING AGENCY NAME(S) AND ADDRESS(ES) 

National Aeronautics and Space Administration 

Washington. DC 20546-0001 

and 

U.S. Army Research Laboratory 
Adelphi, Maryland 20783-1145 


11. SUPPLEMENTARY NOTES 


WU-581-30-13-00 
1L162211A47A 


8. PERFORMING ORGANIZATION 
REPORT NUMBER 


E-12347 


10. SPONSORING/MONITORING 
AGENCY REPORT NUMBER 

NASA TM—2000-210227 
ARL-TR-2246 
DETC2000/PTG-14386 


Prepared for the 8th International Power Transmission and Gearing Conference sponsored by the .American Society of Mechanical 
Engineers, Baltimore. Maryland. September 10-13. 2000. David G. Lewicki and Robert F. Handschuh. U.S. Army Research Labora¬ 
tory; NASA Glenn Research Center: Lisa E. Spievak. Paul A. Wawrzynek, and Anthony R. Ingraffea. Cornell Fracture Group. Cornell 
University; Ithaca, New York 14853-4)001. Responsible person. David G. Lewicki. organization code 5950. (216) 433-3970. 


12a. DISTRIBUTION/AVAILABILITY STATEMENT 

Unclassified - Unlimited 
Subject Category: 37 


12b. DISTRIBUTION CODE 


Distribution: Nonstandard 


This publication is available from the NASA Center for AeroSpace Information. (301) 621-0390. 


13. ABSTRACT (Maximum 200 words) 

Robust gear designs consider not only crack initiation, but crack propagation trajectories for a fail-safe design. In actual 
gear operation, the magnitude as well as the position of the force changes as the gear rotates through the mesh. A study 
to determine the effect of moving gear tooth load on crack propagation predictions was performed. Two-dimensional 
analysis of an involute spur gear and three-dimensional analysis of a spiral-bevel pinion gear using the finite element 
method and boundary element method were studied and compared to experiments. A modified theory for predicting gear 
crack propagation paths based on the criteria of Erdogan and Sih was investigated Crack simulation based on calculated 
stress intensity factors and mixed mode crack angle prediction techniques using a simple static analy sis in which the tooth 
load was located at the highest point of single tooth contact was validated. For three-dimensional analysis, however, the 
analysis was valid only as long as the crack did not approach the contact region on the tooth. 


14. SUBJECT TERMS 


Gears: Spiral-bevel gears; Crack propagation: Fracture mechanics; 
Finite element method; Boundary element method 


17. SECURITY CLASSIFICATION 
OF REPORT 

Unclassified 


18. SECURITY CLASSIFICATION 
OF THIS PAGE 

Unclassified 


19. SECURITY CLASSIFICATION 
OF ABSTRACT 

Unclassified 


15. NUMBER OF PAGES 

15 


16. PRICE CODE 

A03 


20. LIMITATION OF ABSTRACT 


Standard Form 298 (Rev. 2-89) 

Prescribed by ANSI Std. Z39-18 
298-102 


















